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Surface granular flow, comprising granular material flowing on the surface of a heap of 
the same material, occurs in several industrial and natural systems. The rheology of such 
flow is investigated by means of measurements of velocity and number density profiles 
in a quasi- two-dimensional rotating cylinder, half- filled with a model granular material - 
mono-size spherical stainless steel particles. The measurements are made at the center of 
the cylinder, where the flow is fully-developed, using streakline photography and image 
analysis . The stress profile is computed from the number density profile using a force 
balance which takes into account wall friction. Mean velocity and root mean square 
(r.m.s.) velocity profiles are reported for different particle sizes and cylinder rotation 
speeds. The profiles for the mean velocity superimpose when distance is scaled by the 
particle diameter (d) and the velocity by a characteristic shear rate (7c = [gsin(/3 m — 
/3 s )/dcos/3 s ] 1 / 2 ) and the particle diameter where (3 m is the static angle of friction and (3 S 
is the sliding angle of friction. The scaling is also found to work for the r.m.s. velocity 
profiles. The mean velocity is found to decay exponentially with depth in the bed with a 
decay length of A = I .Id. The r.m.s. velocity shows similar behaviour but with A = 1.7d. 
The r.m.s. velocity profile shows two regimes: near the free surface the r.m.s. velocity is 
nearly constant and below a transition point it decays linearly with depth. The shear rate, 
obtained by numerical differentiation of the velocity profile, is not constant anywhere 
in the layer and shows a maximum which occurs at the same depth as the transition 
in the r.m.s. velocity profile. Above the transition point the velocity distributions are 
Gaussian and below the transition point the velocity distributions gradually approach a 
Poisson distribution. The shear stress increases roughly linearly with depth. However, the 
variation of the apparent viscosity (77) with r.m.s. velocity (it) shows a relatively sharp 
transition at the shear rate maximum, and in the region below this point the apparent 
viscosity varies as r\ ~ u~ 1,5 . The measurements indicate that the flow comprises two 
layers: an upper low viscosity layer with a nearly constant r.m.s. velocity and a lower 
layer of increasing viscosity with a decreasing r.m.s. velocity. The thickness of the upper 
layer depends on the local flow rate and is independent of particle diameter while the 
reverse if found to hold for lower layer thickness. The experimental data is compared to 
predictions of three models for granular flow. 



1. Introduction 

Surface flows of granular materials comprise a shallow layer of particles flowing on 
a fixed bed of the same particles. Such flows are encountered in several industrial op- 
erations as well as in nature, fndu strial examp les appear in the proc essing of mate- 
rial i n systems such as rotary kilns l|Peravlll986(l and tumbling mixers l|Perrv fc Greenl 
Il997l) . in the f ormation of heaps for storage and in the feeding and discharge of hop- 
pers and silos <|Neddermanl If 99^) . Natural examples are the formation of sand dunes 
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and transport of sediment in rivers (Bounhourc et al. 2002). 



Surface granular flows have one unique feature that differentiates them from other 
granular shear flows. The 'boundary' between the flowing particles and the fixed bed 
is determined by the flow itself. A consequence is that interchange of particles between 
the flowing layer and the bed is possible, and the local surface angle and the local layer 
thickness vary with the local flow. Although it is reasonable to picture the flow in terms 
of two distinct regions (flowing layer and fixed bed) separated by a sharp interface, in 
fact the boundary is diffuse and the flow decays to zero exponentially with depth. Thus 
a homogeneous and responsive boundary exists at the base of the flow which is quite 
different in qualitative terms from a flow on an inclined rough surface. Several studies 
have been carried out to understand the physics of surface granular flows and a review 
of previous work is given below. 

Rotating cylinders and heaps have been the primary systems used for the experi- 
mental study of surface granular flows. In both systems, the major component of the 
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showed the decay of the velocity to be exponential with depth and 



12002 at lLongo & : Lambertil l2002t ITaberlet , et all 120031: lOrpe fc Khakharfl2004|) . The re- 

gion of exponential decay is relatively small and there has been some focus on the 
velocity gradient corresponding to the linear part. Some studies show that the shear 
rate is nearly constant and independen t of the local mass flow rate, particularly for low 
mass flow rates and near 2d systems l|Raichenbach et oil Il995t [Khakha r et al\ 12001 at 
nam v et all 12002 al). In other cases the shear rate is found to increase with flow rate 
(JOrpe fc Khakhaill200lt I.Tain et all 120021) . The bulk density is nearly constant across 
the layer and approaches close packing ( Raichcnbac h et al 1 11 99,4 Ifiouamv et qj]l2002al 
I.Iain et aZll2002HOroe fc Khakhaill2004l) . indicating the highly dissipative nature of the 
flow. 

The surface angle (0) is foun d to increase with mass flow rate in the layer l|Khakhar et all 
12001 at lOrpe fc Khakhar 20011), however, in these and similar studies the measurements 
are made near a wall. Recently ITaberlet et all (|2003ll have shown a systematic effect of 
the wall on the surface angle. They found that the measured increase in angle with flow 
rate in heap flow experiments could be accounted for by considering the effect of wall 
friction in a force balance. The analysis yields tan/3 = fii + [i w 5/b, where S is the layer 
thickness, b is the gap width, \Xi is the effective coefficient of friction and fx w the coefficient 
of friction beteween the particles and wall. Thus for sufficiently wide g aps (S/b >> 1) a 
constant surface angle should be obtained. This was indeed observed bv l.Iop et all I^OOfij) 
in a systematic experimental study of the effect of gap width on heap flow. They found 
that the gap width between side walls also affects the mean velocity and the depth of 
the flowing layer. As the gap width is increased the surface angle decreases, the mean 
velocity decreases and the depth of the flowing layer increases. The rate of increase of 
surface angle with flow rate reduces with gap width, and at the highest gap width con- 
sidered (570 particle diameters) the surface angle is nearly independent of the flow rate. 
Down-channel velocity profiles measured at the free surface indicate that the range of the 
"wall effect" to extends to a greater distance from each wall with increasing gap width, 
even when the lengths are scaled by the gap width. This, however, may be a reflection o f 
the fact the flow rate varies significantly across the gap. ICourrech du Pont et all 1 2005fl 
measured the instantaneous velocity profiles during an avalanche in a rotating cylinder. 
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The profiles obtained at different times at the centre of the cylinder are self-similar and 
have an exponential dep endence on the d epth. The linear part of the profile, found for 
steady flows, is not seen. Ijain et all l|2002f) have reported measurements of the granular 
temperature, which is an important input in some models for the flow. The qualitative 
features of the flow remain unchanged eve n when the particles are completely immersed 
in a liquid l|Jain. Ottino fc Lueptowll2004|) . 

The simplest models for surface flow comprise a depth averaged mass balance equation 
in the flowing layer with a phenomenological equation for the flu x from the bed to the 
layer assumed to be determined by the local angle of the surfa ce ijBouchaud et all\l99& 
lMehtall99llBoutreux fc de Gennesll996tlBoutreux et qjl998ft. More detailed depth av- 



eraged continuum models als o yield similar equations (Khak har et ah 1997; Douadv^^zZ 
1999t i Khakhar et all 12001 al) and agree well with experimental results IjKhakhar et al 
2001 al 



2001ala ). The main difference in these models is the stress constitutive equation used. 



One approach is to assume the shear ra te to be a known constant which eliminates the 
need for the stress constitutive equation l|Douadv e^"aDll999t lMakselll999j) . rRaichenbachl 
(2003) gave a physical justification for a constant shear rate based on an assumption 
that a single collision entirely dissipates the momentum of both colliding particles. This 
yields an estimate for the shear rate as 7 ~ (gsin/^/d) 1 / 2 , where (3 is the local sur- 
face angle of flow, d is the particle diameter and g is acceleration due to gravity. Such 
an approach gives goo d predictions of experimental results for rotating cylinder flow at 
low rotational speeds llMakselll999l) . A ph enomenological equation for the stress based 
on frictional and colli sional contributions llKhakhar et o./Jll 997t lo7r>e fc Kha,kna,rll200ll 
Kha khar et aiJEoOlol) gives a somewhat different expression for the shear rate 
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g sin(fi m - /3s) " 1 l '~ 
cd cos p s 



(1.1) 



where (3 m is the angle of static friction, (3 S is the static angle of repose (angle of dynamic 
friction) and c is a constant. From the analysis it can be shown that the free surface angle 
for a steady fully developed flow with no interchange is equal to f3 m . The continuum 
model based on the phenomenological stress equatio n gives good prediction s of the layer 
thickness profi le and shear rate for bot h heap flows IjKhakhar et aZ.ll2001a|) and rotating 
cylinder flows I Orpe & Khakhar 200lL The model predicts a nearly constant shear rate 



at low flow rates ^Kh^Mmietal. l20016j) 



More detailed rheological models for the dense gran ular flow, typ ical of surface flows, 
have also been proposed based on different approaches. ISavaee] 1 1998f) developed a kinetic 
theory based on strain-rate fluctuations which were related to the granular temperature. 
The resulting governing equations contain a viscosity term that sho ws a decrease with an 
increase in the granular temperature, similar to that for a liquid. iBocauet et all (|200lj) 
modified granular kinetic theory equations to give a sharper increase in the granular vis- 
cosity with solids volume fraction. The model gave good predictions of the velocity profile 
in a C ouette flow. A two phase solid-liquid model was proposed bv lAranson fc Tsimrind 
( 200^), with the fraction of solid phase present locally determined by an order parameter. 
This model is descri bed in more detail in s ection 13.41 and predictions are compared to 
experimental results. Ijosser and et all ^2004 \ have proposed a continuum model based on 
stress constitutive equations which depend on the local volume fraction. The model gives 
a reasonable qualitative description of surface flows. Dense granular flows have also been 
described by Cosserat type constitutive models in which local couple stress es due to par- 
ticle r otations are taken into account. A Cosserat model was proposed bv lMohan et all 
(2002) and predictions of the model for the velocity profile were shown in good agree- 



4 



Ashish V. Orpe and D. V. Khakhar 



merit with measurements in a vertica l channel flow and dense Couette flow. Recently 
Pouliquen and coworkers l)MiDil l2004t I.Tod. Forterre fc Pouliauenl [2005) proposed that 
the local shear stress in a granular flow may be expressed in terms of a shear rate (7) 
dependent effective friction coefficient of the following form 

/i = /is + V7TI' (L2) 

where / = d^/(p p /P) 1 ^ 2 is the dimensionless shear rate, p v the particle density and P 
the granular pressure. /i2 and Iq are model param eters. The model is found to wor k 
well for flow on a rough inclined plane and heap flow ||Jop. Forterre fc Pouliauenll2005|) . 
We compare our experimental reults to this model as well. 

All the models discussed above are based on a local stress constitutive equation. 
The high density of the flow results in multi-particle contacts and thus there is a dis- 
tinct possibility of non-local transfer of str esses. Several theories have explored such 
ideas. IPouliauen. Forterre fc Le Dizesl l)200l|) considered the stress transfer by fluctua- 
tions in the bed. The fluctuations decay with distance from their point of origin. The 
model gives good predictions of data for dense flows in y erical and inclined channels. 
Mills. Loggia fc Tixierl <|l999h and lBonamv fc Mills! l|2003h have considered dense shear 
flows in which near vertical particle chains are formed due to gravity forces. The model 
gives good predictions of measured velocity profiles. Predictions of this model are also 
compared to our experimental data below. 

The rheology of surface granular flows has not been experimentally explored in detail 
as yet and is the focus of this work. The objective of the present work is to carry out 
an experimental study and analysis of the flow and rheology of a steady, fully-developed 
surface granular flow. A quasi-2-dimensional rotating cylinder system (gap width 5-10 
particle diameters) is used and the measurements are made in the region near the cylinder 
axis where the flow is fully developed and steady. Experiments are carried out over a 
range of particle sizes and rotational speeds, and the data are analyzed by scaling, and 
comparison to existing models. Section gives the details of the experimental system 
and methods used. The results are discussed in section followed by conclusions of the 
work in section 



2. Experimental details 

We first describe the experimental system used for this work. The experimental method 
based on digital photography is discussed next, followed by a description of the technique 
used for analyzing the digital images. 

2.1. System 

Experiments are carried out in quasi- 2-dimensional aluminium cylinders (length b = 1 
or 2 cm) of radius 16 cm. The end walls are made of glass for visualizing the particle 
motion. A computer controlled stepper motor with a sufficiently small step is used to 
rotate the cylinders. Stainless steel (SS) balls of three different sizes with diameters (d) 
1, 2, and 3 mm are used in the experiments. A few experiments are also carried out using 
2 mm brass balls to investigate the effect of material characteristics. Brass is a relatively 
heavier and softer material as compared to SS. The properties of all the particles used are 
listed in tabled All the particles used in the experiments are highly spherical in shape 
(the SS balls are used for manufacturing ball bearings). The surface of all the particles 
is smooth and reflective which facilitates the image analysis method used in this work. 
Experiments are carried out at different rotational speeds (u — 2 — 9 rpm) such that 
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Material Diameter (d) Mass (m) Particle Bulk Young's 

density (p p ) density (pb) Modulus 
mm mg kg/m 3 kg/m 3 GPa 



1 4.2 ± 0.06 4973 

SS 2 33.5 ± 0.18 8000 4899 200 

3 112.5 ± 0.18 4725 

Brass 2 35.0 ± 0.25 8500 5205 120 

Table 1. Geometrical and physical properties of particles. 




Figure 1 . Schematic view of the flow in a rotating cylinder, comprising a surface flowing layer 
of thickness 8{x) and a packed ped of particles rotating with the cylinder. The coordinate system 
employed in the analysis is shown. 

the flow is always in the rolling regime. The typical flow pattern, comprising a surface 
flowing layer of varying thickness (S(x)) and the bed rotating at the angular speed of 
the drum (w), is shown schematically in figure The particles are heavy enough and 
conductive so that static charge effects are negligible. All the experiments are carried out 
with the cylinder half-full of particles (i.e., a fill fraction of 50%). 

2.2. Experimental method 
The measurements in the flowing layer are made near the center of the cylinder (x = 0, in 
figure where the layer thickness is maximum and the flow is nearly unidirectional and 
non- accelerating (velocity is nearly invariant in the cc-direction). The images are taken 
keeping the camera close to the face plate of the cylinder and framing a small region of 
flow. The size of the recorded region is 2560 x 1920 pixels, with one pixel corresponding 
to 0.016 — 0.03 mm depending on the distance of the camera from the cylinder. A point 
source of light is directed nearly parallel to the face plate of the cylinder so as to illuminate 
only the front layer of the flowing particles. Owing to the high sphericity and the polished 
surface of the particles, a well-defined reflection is produced from each stationary particle 
(see figure|2 lower right corner). Each moving particle thus generates a streak of definite 
length depending on its speed and the shutter speed of the camera. Figure |3 shows a 
typical image with streaks of different lengths in the flowing layer. The velocity of each 
particle is determined from the length and the orientation of the streak formed; the 
analysis met hod is expla ined in |2.3I This method of velocity determination is similar to 
that used bv lRaichenbach et all l|1995h . 
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Figure 2. Flow of SS balls of diameter 2 mm for an angular speed (w) of 3 rpm. Streaks 
correspond to the particle displacement during the time the shutter is open At — 1/500 s. x,y 
correspond to the coordinates with the t/-axis along the free surface of the flowing layer, while 
X,Y correspond to the coordinates in the laboratory reference frame. 
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Figure 3. Image analysis of streaks, (a) Original image of the streaks, (b) Thresholded streaks 
with pixels identified in each streak, (c) Fitting of the function 12. 1 1 to one streak. X, Y: co- 
ordinate axis in the laboratory reference frame. X C ,Y C : Position of the streak in this frame of 
reference. 21 and 2w are the length and width of the streak, respectively. 8 is the orientation of 
the streak with the horizontal. 



Images are taken for a range of camera shutter speeds (1/15 — 1/1000 s) so as to account 
for the varying velocity across the flowing layer. About one hundred images are taken for 
each shutter speed with an overall of thousand images combined over different shutter 
speeds. Each experiment is carried out for about 500 revolutions with two images taken 
per revolution. Thus the flow is sampled over a considerable time period (~ 3 h). The 
images are analyzed to determine the velocity and the number density profiles across the 
layer. 

2.3. Image analysis 

The image analysis technique involves identification of the streaks in the front layer 
followed by determination of the length and the orientation of each identified streak. 
Computer codes have been written to manipulate and analyze the digital images. Details 
are given below. 

Each image is first converted to black and white by putting intensity values of pixels 
greater than a threshold to white and the rest to black. The threshold is chosen so 
that the streaks from the front layer are converted to white and the streaks in the back 
layers, which are less illuminated, to black. All the photographs are taken with a black 
background so as to make the thresholding easier. The thresholded image is scanned row 
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by row to identify clusters of contiguous white pixels. Clusters smaller than a specified 
size (75 pixels) are due to reflection of the light from the periphery of the particle surface 
and are eliminated. Clusters larger than this size are taken to be streaks. Figure^ shows 
three magnified streaks in the front layer and a number of faint streaks from the back 
layers. Figure [3p shows the thresholded image with faint streaks eliminated. 

In the next stage of the analysis, the position, length and angle of each streak is 
determined by fitting a parameterized intensity function to the intensity values of the 
streak pixels (and an immediate neighbourhood) in the original image. An intensity 
function of the following form is used 

Ith = Ce-y 2 / Aw2 \x\ < I (2.1a) 

= c e [-s 2 -(2-0 2 ]/4^ 2 | S | > j (2.i 6 ) 

where x = {X - X c ) cos 9 + (Y - A) sin 9, y = {Y - Y c ) cos 9 - (X - X c ) sin 6, 21 is the 
nominal length of the streak, 2w is the nominal width of the streak, 6 is the orientation 
of the streak with the horizontal and C is the maximum intensity value obtained at the 
centroid of the streak. (X c , Y c ) is the position of the centroid of the streak in laboratory 
coordinates (Fig. HJ). 

Initial estimates of the position (X c , Y c ), length (2/), width (2w) and orientation angle 
(9) are first obtained for an identified streak from the thresholded image. The values 
are refined by m inimizing the following function for each streak using Powell's method 
l|Press et a?Jll994|) 

/ = yiC^is/t — I image) 2 1 (2-2) 

where Iimage is the intensity value for a pixel in the original image at that point (fig- 
ure |3]i) . The summation in 1|2.2|1 is carried over all the pixels of the streak together with 
six layers of pixels added to the cluster. Figure shows the computed intensity values 
plotted on the original image using the fitted parameters. 

The velocity magnitude is v = 2L/At, where At is the time interval corresponding to 
the shutter speed set in the camera, and the velocity components are calculated from 
the orientation angle, 9. L = L(l,w) is the distance moved by the particle in the time 
interval At and is obtained as a function of the streak length and width using calibration 
experiments described below. The position corresponding to this velocity vector is taken 
to be the centroid of the streak (X c , Y c ). A similar procedure is applied for every streak 
identified in an image. Only streaks with lengths in the range 20 — 70 pixels are considered 
for analysis. This is because streaks longer than 70 pixels tend to overlap with other 
streaks, while those shorter than 20 pixels result in an inaccurate determination of the 
orientation angle (9). The use of different shutter speeds gives streaks with lengths within 
this range at each depth in the layer. Streaks which are cut off at the edges of an image 
are not considered for velocity calculation. 

The analysis results from all the images are combined and the coordinate positions for 
all the selected streaks are rotated by angle 9 avg . This angle is determined by averaging 
the angle of orientation (9) for all the streaks in the selected region of the flowing layer, 
and represents the orientation of the flowing layer with the horizontal. The components 
of the instantaneous velocity parallel (c x ) and perpendicular (c y ) to the flow direction 
are determined for every streak from its orientation angle 9. The region shown in figure 
is then divided into bins of depth equal to the particle diameter and length 20 mm 
parallel to the flowing layer. The profiles of the averaged parallel (v x ) and perpendicular 
(v v ) components of velocity are then generated by averaging the instantaneous velocities 
(c x , c y ) over all streaks in each bin. 
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In a typical measured velocity profile, the average a> velocity (v x ) decreases and then 
starts to increase with depth in the flowing layer from the free surface. This is because 
velocity decreases and eventually becomes negative (due to bed rotation) with decreasing 
y. However, our technique cannot give the direction of the velocity vector. Hence we 
assume all velocities below the minimum velocity point to be negative. 

The root mean square (r.m.s.) velocity and its components in x and y direction are 
calculated from the instantaneous velocities using following equations 



where (©) represents an average over all streaks in a bin corresponding to a y-coordinatc. 
The number density profiles are obtained from a separate set of images at a high shutter 
speed (1/2000 s). Images taken at such high shutter speeds ensure that all streaks are 
short (< 30 pixels) and do not overlap so that all streaks in an image can be used for 
analysis. This is required to get an accurate number count. Images are analyzed in the 
same manner as mentioned above, but only the centroid positions are noted. The number 
density for any y is determined by counting the number of particles in each bin. A source 
of inaccuracy is that particles from back layers may also be counted which leads to an 
overestimate of the number densities. This is particularly the case for the low number 
density region near the free surface. Further, the measurements represent the densities 
near the wall which are different from those in the bulk. 

Experiments are carried out to calibrate the analysis technique. A single particle is 
glued to the face plate of the cylinder from the inner side at a position close to the 
periphery of the cylinder. The radial position of the particle is accurately determined 
from which the velocity of the particle can be obtained since the angular speed of the 
cylinder is known. The cylinder is rotated at different angular speeds and images of 
this particle are taken at different camera shutter speeds (1/30 s - 1/500 s). The size 
of the recorded region and the lighting conditions are the same as those used for the 
measurements in the flowing layer. Further, the rotational speed of the cylinder and 
the shutter speeds are adjusted so as to generate streaks of three different lengths. The 
maximum and the minimum lengths of the streaks span the entire range of streak lengths 
considered for velocity measurements in the flowing layer. The experiments are carried 
out for all the particles used in the experiments. The length (21), width (2w), position 
(X C ,Y C ) and orientation (6) of this single streak are determined in each case using the 
same analysis procedure as described earlier. The actual distance traveled by the particle 
at radial distance r (L = ujrAt) is known from the linear velocity of the particle (tor) 
and the shutter speed used (At). A linear relation of the following form is assumed to 
hold between the actual length (L) and measured length (I) and width (w) of each streak 



The constants a\ and ai are determined by a least-squares fit to the experimental data 
combined for particles of different sizes and materials. The fit yields a\ — 2.05 and 
a.2 = 1.86 with the standard error of about 1.1 pixels. Figure^ shows a comparison of 
the measured and fitted streaks. The experimental data are well correlated by (|2.4(l and 
indicate that the actual distance traveled by a particle depends on both the length and 
width parameters of the intensity function. The error in calculating L is less than 3% 




(2.3a) 



(2.36) 



(2.3c) 



L = ail + a2iv. 



(2.4) 
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Figure 4. Calibration results for determining streak lengths. Variation of actual streak length 
(L) with a linear sum of length (I) and width (w) of a streak measured by image analysis. 
a± — 2.05 and 02 = 1.86 are the constants obtained by linear regression of the experimental 
data. Dotted lines represent the 3% deviation in calculation of L. 
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Figure 5. Mean velocity profiles measured at three locations along the layer as indicated in 
the legend, for 2 mm particles with the cylinder rotated at 3 rpm. 



(figure |2J and represents the typical error in measurements of velocities in the flowing 
layer. 

The dynamic angle of repose (fl m ) is taken to be the angle by which the data is rotated 
{Oavg) to obtain v y — 0. The static angle of repose (f3 s ) is taken to be the angle of the 
stationary free surface when the rotation of the cylinder, operating under steady flow 
conditions, is stopped. The values of /3 m and j3 s (averaged over six images) are reported 
in table |21 for all the cases studied. We also checked that the flow was fully developed 
by measuring the velocity profile at two additional locations: one 10 mm upstream and 
one 10 mm downstream of the layer midpoint [x — 0). Figure [3] shows the profiles of the 
mean velocity at the three locations for one case (2 mm particles in the cylinder rotated 
at 3 rpm). The profiles are identical within experimental error, verifying the the flow is 
fully-developed. 
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Material d u> /3 m f3 3 q 

(mm) (rpm) (deg.) (deg.) (cm 2 /s) 





1 


2 


25.6 


16.0 


26.8 




1 


3 


26.6 


16.0 


40.2 


SS 


2 


3 


29.1 


20.1 






2 


6 


35.7 


20.1 


80.4 




2 


9 


42.8 


20.1 


120.6 




3 


3 


26.3 


18.4 




Brass 


2 


3 


27.8 


18.3 





Table 2. Experimentally determined values of the angles /3 m and f3 3 , and calculated 
volumetric flow rate per unit width (q) for the different cases studied. 



3. Results and discussion 

We first present the experimental results obtained for varying system parameters. 
Scaling and the analysis of the experimental profiles are presented next. The rheology of 
the flow is discussed in 13.31 Finally, we compa r e our e xperimental velocity profiles with 
predictions of the models of iBonamv fc Mills! (12003 ). lAranson fc Tsimrind l)2002|) and 
Pouliqucn and cowor kers |(,Tod. Forterre fc PouliaiierfeOO^ lMiDill2004lL 

3.1. Base data 

The primary parameters varied in the experiments are the cylinder rotational speed (w) 
and the particle diameter (d), and the results presented below are in terms of these 
parameters. There is, however, a direct relation between the cylinder rotational speed 
and the local volumetric flow rate per unit gap width (q) . A mass balance for a half- full 
cylinder, assuming that the bulk density in the bed and the layer are the same, yields 
q = ojR 2 /2, and calculated values are given in tabled Increasing the rotational speed is 
thus equivalent to increasing the local flow rate. 

The variation of the x component of the mean velocity (v x ) with depth (y) in the 
layer is shown in figure for different rotational speeds, particle sizes and materials 
studied. The error bars denote the standard deviation over ten different data sets. All 
the measured quantities are plotted from y = at the free surface of the flow down to 
about 3-4 particle diameters below the point corresponding to the minimum absolute 
velocity. The free surface is taken to be the point at which the measured number density 
is close to zero (n < 10 -3 mm -2 ). The profiles are linear over most of the layer depth 
for all the cases studied, with an exponential decay near the base of the flowing layer. 
The negative velocities correspond to the particles in the rotating bed which move in a 
direction opposite to that of the flowing layer. Further, the profiles tend to flatten out 
near the free surface. This is a low density region comprising saltating particles. The 
ballistic trajectories in this region are clearly seen in figure Q which is an image of the 
flow taken at a low shutter speed. 

The magnitude of the maximum velocity and the depth of the layer increase with the 
rotational speed for a fixed particle size (figureUJi). The length of central linear region of 
the velocity profile increases with increasing rotational speed, while the length of the ex- 
ponential region near the bottom of the layer is nearly constant for all the cases. The shear 
rate (7 = dv x /dy) corresponding to the linear portion increases with rotational speed. 
The profiles for SS and brass particles are close to each other indicating that particle 
material parameters like coefficient of restitution and particle roughness have a minor 
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Figure 6. Variation of a:-component of mean velocity (v x ) with depth (y) in the layer, (a) Thick 
solid line denote experimental data for 2 mm brass particles rotated at 3 rpm. Other lines denote 
data for 2 mm SS balls at different rotational speeds. Thick dotted line shows a straight line 
fitted to the linear portion of the 9 rpm profile, (b) Thick solid line denote data for 1 mm SS 
balls rotated at 2 rpm. Other lines denote data for SS balls rotated at 3 rpm. Error bars denote 
standard deviation over 10 sets. 



Figure 7. Typical time lapse image showing ballistic trajectories of particles in the low 
density region (above the dashed line). 

effect on t he flow. Th i s is in agreement with our earlier work llOrne fc Khakharll200lh 
and that of lJain et jU\ {2002) in rotating cylinders. The maximum velocity decreases and 
the depth of flowing layer increases with increasing particle diameter (figure ^p) . The 
shear rate, thus, decreases with increasing particle size. 

Qualitatively similar velocity profiles have been reported in the literature for different 
experimental systems. A linear and an exponential region of ve locity profile was ob- 
served in rotating cylinders for various rotational spe eds studied (Bo namv et a/.ll20 02a; 
lLongo fe Lambertll2002t IZanuttigh fe Lambertill200^) . T he trends obtaine d in figure 03 
for varying particle size are similar to those reported bv I.Tain et all l(2002f ). The system 
size used in all these cases is comparable to our experimental system. The velocity pro- 
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Figure 8. Variation of y-component of mean velocity (v y ) with depth (y) in the flowing layer, 
(a) Open symbols denote experimental data for 2 mm SS balls. Filled symbols denote data for 
2 mm brass particles rotated at 3 rpm. (b) Open symbols denote experimental data for SS balls 
for a rotational speed of 3 rpm. Filled symbols denote data for 1 mm SS balls rotated at 2 rpm. 
Error bars denote standard deviation over 10 sets. 
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Figure 9. Variation of area fraction (<f>) with depth (y) in the flowing layer, (a) Open symbols 
denote experimental data for 2 mm SS balls. Filled symbols denote data for 2 mm brass particles 
rotated at 3 rpm. (b) Open symbols denote experimental data for SS balls for a rotational speed 
of 3 rpm. Filled symbols denote data for 1 mm SS balls rotated at 2 rpm. Error bars denote 
standard deviation over 10 sets. 



file i n the lower ha lf of th e layer resembles the profiles obtained bv lMueth et all l)200fj|) 
and iBoc quet et al\ (|200l f) for Couette flows. I n some of the studies mentioned above 
ijRaichenbach et a/lll995HBonamv et alJl2002a|h the flattened upper low density region 
was not obtained. 

The magnitude of the y-component of the mean velocity (v y ) is close to zero throughout 
the depth of the layer for all the cases (figure [SJ). This confirms that the flow is nearly 
unidirectional (invariant in the x direction) in the region of measurement. The profiles 
of the projected area fraction of the particles <p = mrd 2 /4A, where n is the number of 
particles counted per bin and A is the area of the bin are shown in figure The area 
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Figure 10. Variation of the components of r.m.s. velocity with depth in the flowing layer. Open 
symbols denote experimental data for 2 mm SS balls. Filled symbols denote experimental data 
for 2 mm brass particles rotated at 3 rpm. Error bars denote standard deviation over 10 sets. 
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Figure 11. Variation of the components of r.m.s. velocity with depth in the flowing layer. Open 
symbols denote data for SS balls rotated at 3 rpm. Filled symbols denote data for 1 mm SS 
balls rotated at 2 rpm. Error bars denote standard deviation over 10 sets. 



fraction is nearly constant over the entire layer for all the particles with a rapid drop 
near the free surface. No appreciable change in the magnitude is observed for different 
rotational speeds and materials. The area fractions are significantly smaller than the 
hexagonal close packed limit (4> c = 0.866). Nearly constant s urface densities acros s the 
flowing layer wer e also observed bv lRaichenbach et al\ l)l995|) . lBonamv et all l|200lj) and 
IJain et all <|2002() . 

The components of the r.m.s. velocities {u x ,u v ) across the flowing layer are shown in 
figures ITUI and ITTI In both cases, the r.m.s. velocity profiles are nearly linear over a large 
portion of the layer with a flattening near the free surface and an exponential decay 
near the base of the flowing layer. The sizes of the flattened region as well as the region 
of exponential decay are larger for the r.m.s. velocity profiles than the corresponding 
regions for the mean velocity profiles. The r.m.s. velocities are smaller by nearly an order 
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Figure 12. Variation of the ratio of the r.m.s. velocity components across the flowing layer, 
(a) Open symbols denote data for 2 mm SS balls. Filled symbols denote data for 2 mm brass 
particles, (b) Open symbols denote data for SS balls rotated at 3 rpm. Filled symbols denote 
data for 1 mm SS balls rotated at 2 rpm. 



of magnitude as compared to the mean velocities. The relative error is thus larger in 
this case as compared to the case of mean velocities, as indicated by the error bars. The 
r.m.s. velocities increase with rotational speed for a fixed particle size. The slope of the 
linear region in each case is nearly independent of the rotational speed, but increases 
with particle diameter (i.e., du x /dy,du y /dy decrease). The r.m.s. velocities for the brass 
balls are quite close to those for steel which again indicates that the flow behaviour is 
nearly independent of the material properties. The anisotropy in the r.m.s. velocities is 
clearly evident from the graphs and the r.m.s. velocity component in the flow direction 
is larger by about 50% than that normal to the flow directio n. A linear and a flattened 
region in the r.m.s. velocity profiles was also obtained by l.lain et all l)2002t) in their 
experiments in rotating cylinders but they found the r.m.s. velocities to be about l/3rd 
of the mean velocity near the free surface. These values are significantly larger than the 
r.m.s. velocities obtained here. 

The degree of anisotropy in the r.m.s. velocities across the flowing layer for all the 
cases studied is shown in figure IP21 The ratio u x /u y (« 1.5) is nearly constant across 
the flowing layer and independent of rotational speed, while it decreases slightly with an 
increase in the particle diamete r . The data for brass and SS balls is nearly the same. 
lAzanza. Chevoir fc Moucherontl l)l999|) found the ratio of the r.m.s. velocity compo- 
nents to be 1 . 45 for 2-dimensional chute flows, while a ratio of 1.3 was obtained by 
iBocauet et all ()200lf) for Couet t e flow s. Anisotropy in the r.m.s. velocity components 
was also observed bv l.lain et all l)2002|) for the experiments in rotating cylinders. 

We obtain the stress profiles in the flowing layer from a stress balance taking into 
account the contribution of wall friction. Assuming the flow to be steady and fully de- 
veloped, the Cauchy equations averaged across the gap width (b) are 

= + I T zx \l + (p)g sin f3 m (3.1a) 

ay 

where (•) = (1/6) J Q -dz denotes an average across the gap. The stress at the side walls 
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due to friction is 

T zx \o = — T Z x\b = TzzHw, (3-2) 

where p w is the coefficient of wall friction. The coefficient of wall friction is assumed be a 
constant, which may overestimate the actual value. Assuming further that the variation 
of stresses across the gap is small, that the normal stresses are isotropic (t zz = T yy ) and 
that the stresses vanish at the free surface we obtain 

,o 

Tyy = pb cos m / (4>/4>b)dy (3.3a) 
J —y 

(<p/(f) b )dy - 2p w /b / Tyydy. (3.36) 
-v J-v 

In the above equations 4>b is taken to be the constant value near the base of the flowing 
layer (figure|5J). The bulk density in the rotating bed is obtained as p b = l/2(rnb/irR 2 b), 
from the measured mass of particles in the half-filled cylinder, mj, and results are given 
in table ^ The effective coefficient of friction for the particles at the base of the layer is 
obtained from l|3.3|l as 

Pi = — = tan/3 m - fhv^r, (3.4) 

Tyy 

where 8 e = [2 T yydy\/[p b f_ y {4>/4>b)dy\ is the effective layer thickness and y s is the 
position of the base of the layer, which we take to be the point a t which the shear rate 
extrapolates to zero. Equation l|3.4|l is identical to one obtained bv lTaberlet et aj\ l)2003[) 
if 5 e is replaced by the measured layer thickness (5). Here we obtain S e by numerical 
intergration of the experimental profiles to avoid the ambiguity in defining the location 
of the free surface. 

Equation l|3.4|) allows the estimation of the wall friction coefficient from experimental 
data. FigureElshows a plot of the measured values of tan (3 m versus S e /b for all the cases 
studied. All the data at 3 rpm (open symols), corresponding to a fixed flow rate, fall on 
a straight line which yields a wall friction coefficient of p w = 0.09, which is a reasonable 
value for steel beads on a smooth surface. The data for higher rotational speeds (filled 
symbols), which correspond to significantly higher flow rates, deviate from the fitted 
straight line. This seems to imply t hat the effective coe fficient of friction (pi) may be 
velocity dependent, as suggested bv lDouadv et al\ (^.999). In the data reported here we 
correct the shear stress using a wall friction coefficient of p w = 0.09 in all cases. 

The profiles for the shear stress are shown in figure 1711 The stresses are linear across 
most of the layer with a slight curvature near the free surface and a deviation from 
linearity deeper in the layer due to wall friction. The linear variation results from the 
nearly constant area fraction (</>) over most of the layer. The shear stress increases with 
the angular speed (figure Et), while the normal stress shows the opposite trend. This is 
because (3 m increases with angular speed. The shear stress profile is nearly independent 
of the particle size (d) and the material type (figure El 3 )- 

The compone nts of the streaming stress contribution to the total stress tensor are 
JCamnbelJll99(ih 

T yx(v) = p(c^y~~ C^Cyl (3.5a) 
T xx(y) =p(c^-c^c r ) (3.56) 
Tyyiv) = p(cyCy~-Cy~~5i) (3-5c) 
where c x and c y are the components of the instantaneous velocities for all the relevant 
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Figure 13. Variation of coefficient of static friction (tan/3 m ) with scaled effective layer thickness 
(Se). Open symbols: data for 3 rpm. Filled symbols: data for 6 and 9 rpm. Line is least squares 
fit to 3 rpm data. 
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Figure 14. Shear stress (r yx ) profiles in the flowing layer, (a) Open symbols denote data for 2 
mm SS balls. Filled symbols denote data for 2 mm brass particles rotated at 3 rpm. (b) Open 
symbols denote data for SS balls rotated at 3 rpm. Filled symbols denote data for 1 mm SS 
balls rotated at 2 rpm. 



particles for every y-coordinate and p is the nearly constant layer density. The profiles 
of the streaming component of the shear stress normalized by the total shear stress are 
shown in figure iFol Very similar profiles are obtained for the normal stresses. Streaming 
stresses are a significant component of the total stress near the free surface, and are more 
than 20% of the total stress in the region corresponding roughly to the flattened portion 
in r.m.s. velocity profiles. Deeper in the bed the streaming stresses become negligibly 
small. This behaviour is consistent with the observed particle trajectories (figured) and 
the density profile (figure El . 
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Figure 15. Variation of the ratio of streaming stress (Ty X ) to the total shear stress (r yx ) with 
depth (y) in the flowing layer, (a) Open symbols denote data for 2 mm SS balls. Filled circles 
denote data for 2 mm brass particles rotated at 3 rpm. (b) Open symbols denote data for SS 
balls rotated at 3 rpm. Filled triangles denote data for 1 mm SS balls rotated at 2 rpm. 



3.2. Scaling and analysis 

We rescale the velocity profiles using a characteristic velocity vc = Jed, with the char- 
acteristic shear rate given jc = [gsm((p m — /3)/dcos/3 s ] 1,/2 based on l|l.ljl . and use the 
particle diameter (d) as the characteristic length. After rescaling, the dimensionless pro- 
files are shifted in the y direction so as to superimpose the profiles. The resulting profiles 
are shown in figure HBk . All the profiles for different rotational speeds, particle sizes and 
material types collapse to a single curve. The linear regions of all the profiles, in particu- 
lar, superimpose very well. The scaling thus describes the variation of the mean velocity 
profiles with particle diameter and local flow rate quite well. Figure ITfib shows the same 
data but reseated b y the characteristic velocity based on the shear rate proposed by 
iRaichenbachl l)2003j) . ja = (g sin j3 m / d) 1 ' ' 2 . The scaled profiles superimpose equally well 
in this case, after appropriately translating the profiles in the y direction. From the data 
available it is difficult to distinguish between the two relations for the shear rate. 

The r.m.s. velocity profiles scaled using the characteristic shear rate jc and shifted 
by the same distance as that for the mean velocity profiles are shown in figure ITTk . The 
scaled profiles collapse to a single curve over most of the layer except for the region near 
the free surface, though the relative scatter in this case is greater as compared to the 
mean velocity profile. 

We define an average shear rate for a velocity profile, j a , as the velocity gradient 
obtained by fitting a straight line to the linear portion of the velocity profile. An example 
of such a fit is shown in figure We note that the magnitude of the shear rate varies 
depending on the number of points included in the fitting. This variation is, however, 
less than 10% when the number of points that are assumed to lie in the "linear" portion 
are varied. Figure 1181 shows a comparison between the measured shear rates and the 
characteristic shear rate jc- Although there is scatter in the data, there is reasonable 
agreement between the model predictions and expe riment. A least squares fi t gives c = 
0.56 which is different from the value estimated bv lOrpe fc Khakharl 1(2001) (c rs 1.5). 
The latter value was based on shear rates calculated assuming a linear velocity profile 
across the entire flowing layer. The velocity profiles are linear only in the central region of 
the layer and the shear rates are lower near the free surface and the bed. The shear rates 
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Figure 16. Variation of the rescaled mean velocity with rescaled depth in the flowing layer. 
Filled symbols denote data for brass particles rotated at 3 rpm. The data for all sets are shifted 
after scaling so as to superimpos e the profiles, (a) Sc aling based on shear rate given by 11.11 . 
(b) Scaling based on shear rate of Raichenbach (2003) 
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Figure 17. Scaled profiles of the components of r.m.s. velocity in the flowing layer. Filled 
symbols denote data for brass particles rotated at 3 rpm. 



calculated bv lOrpe fc Khakharl fcOOll) are, thus, lower than the values corresponding to 
the central linear part, which results in a higher value for c. 

The average shear rate (7 a ) is also proportional to the shear rate proposed bv lRa i chenbachl 
(2003) (tr), as shown in figure [THb. However, a linear fit gives a significant negative in- 
tercept, with 7 a = 1.26(tr — 16.9). This comparison indicates that -jc ma Y give a better 
description of the flow than tr. 

A semi-log plot of the velocity profile indicates that the exponential region of the 
velocity profile spans 3-4 particle diameters ( figure m?|l . We fit straight lines to this data 
to obtain the characteristic decay length, A, defined by 



(3.6) 
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Figu re 18. (a) Comparison of the measured average shear rate (%) to the shear rate predicted 
from iKhakhar etaH (12001 oft (-y e ). (b) Comp arison of the measured average shear rate ( A / a ) to 
the shear rate predicted from iRaichenbachl l|2003f) (tr). Symbols show the experimental data 
for all the cases studied. Solid lines show linear fits to the data. 

where vq is a constant of proportionality. The variation of the characteristic length A 
with the particle diameter (d) for the different cases studied is shown in figure The 
values of A depend on the particle size (d) , but are nearly independent of rotational speed 
and material type. A linear fit to the data yields X/d = 1.1 ± 0.2. This result indicates a 
relatively sharp decay of the velocity with depth in the bed since the characteristic length 
is just a little more than one particle diameter. A l inear variation o f the characteristic 
decay length with particle diameter was obtained bv lKomatsu et a,l\ l)200ll) for exponen- 
tia l region in heap flow . A significantly larger characteristic length (10d) was obtained 
bv lMohan et all |l997) for dense flows in a vertical channel. 

An exponential decay is also found for the portion of the r.m.s. velocities near the 
base of the layer. A plot of the characteristic length with the particle diameter obtained 
from fits to the r.m.s. velocity (u) is shown in figure l2TJb. A linear fit to the data yields 
X/d = 1.7±0.3. Thus the r.m.s. velocity decays slightly slower with depth than the mean 
velocity. 

3.3. Rheology 

The variation of the measured quantities with depth in the layer for one case are shown 
side-by-side in figure l^Tl to give an overall qualitative view of the rheology. The shear rate 
profile, obtained by numerical differentiation of the mean velocity profile (7(2/) = dv x /dy), 
is shown in figure For all the cases studied, the shear rate increases to a maximum 
with increasing depth and then decreases as the fixed bed is approached. We note that a 
relatively smooth curve is obtained as a result of the averaging with the bin size equal to 
the pa rticle diameter. An oscillatin g shear rate pr ofile is obtained, as reported in previous 
works l|Mueth et alJl200fi lMuethll2003l Trim et a/jboolh . if a smaller bin size is used. The 
fluctuations in the shear rate profile arise from the magnification of errors as a result 
of taking a derivative of the velocity profile data. The fluctuations are larger for higher 
rotational speeds and smaller particles. The results indicate that the velocity profile is 
not linear anywhere in the flowing layer, although the region around the maximum may 
be approximated to be linear. The shear rate decreases to a small value near the free 
surface. 
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Figure 19. Mean velocity (u^) profiles across the flowing layer plotted on a semi-log scale. Solid 
lines show exponential fit to the linear portion of the profiles, (a) Open symbols denote data for 
2 mm SS balls. Filled circles denote data for 2 mm brass particles rotated at 3 rpm. (b) Open 
symbols denote data for SS balls rotated at 3 rpm. Filled triangles denote data for 1 mm SS 
balls rotated at 2 rpm. 
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Figure 20. Variation of the characteristic length (y*) in the exponential fit with the particle 
diameter (d) for (a) mean velocity profile (b) r.m.s. velocity profile. 



The graphs in figure |^ indicate that on the basis of rheological behaviour, the layer 
can be considered to comprise two regions with different flow characteristics. Above a 
critical depth, y c , corresponding to the depth of the shear rate maximum, the shear stress 
increases with shear rate. However, in the region below y c , the trend is reversed and the 
shear stress increases while the shear rate decreases. There also appears to be a transition 
in the r.m.s. velocity at y c : above y c the profile is nearly fiat, whereas it decreases with 
depth below y c . Similar results are obtained for all cases studied. 

Figure 1221 shows a comparison of profiles for 2 mm particles and different cylinder 
rotational speeds. All graphs for a particular rotational speed are shifted along the y axis 
so that the area fraction profiles superimpose on each other. This ensures that the free 
surface in all cases is located at the same y position so that depths may be compared. 
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Figure 21. Profiles of the mean velocity (v x ), r.m.s. velocity (u), area fraction (<j>), shear rate 
(7) and shear stress (t vx ) in the flowing layer for 2 mm SS balls rotated at 6 rpm. Dashed line 
denotes the transition depth (y c ). 




Figure 22. Profiles of the mean velocity (v x ), r.m.s. velocity (u), area fraction (<f>), shear rate 
(7) and shear stress {r yx ) in the flowing layer for 2 mm SS balls at different rotational speeds. 
Horizontal lines denote the transition depth for each rotational speed. 



The transition in the flow is seen in all the cases and the data indicate that y c increases 
with rotational speed. Figure [231 shows a comparison of profiles for different diameter 
particles rotated at 3 rpm. In this case the depth is rescaled by the particle diameter 
before shifting the graphs to superimpose the area fraction profiles. The results indicate 
that the scaled depth decreases with increasing particle diameter. Finally, figure l2~H shows 
data for 3 mm particles for two different gap widths (6) . Although the mean velocity and 
r.m.s. velocity profiles are slightly different for the two gaps, the transition depth does 
not depend on the gap width. 

Evidence for transition in the flow is also seen in the velocity distributions. Figure 12*51 
shows the velocity distributions at different locations in the layer for 2 mm particles 
rotated at 3 rpm. The distributions vary significantly with depth in the layer. The c y dis- 
tribution is Gaussian for all points above the transition point (y > y c ) ffigurel2"5b). Below 
the transition point (y < y c ) it gradually evolves to a Poisson distribution with increasing 





Figure 24. Profiles of the mean velocity (v x ), r.m.s. velocity (u), area fraction (<j>), shear rate 
(7) and shear stress {r yx ) across the flowing layer for 3 mm SS balls at 3 rpm for different gap 
widths (6). The dashed line denotes the transition depth. 



depth in the layer ffigure ESt O. The c x distribution is Gaussian above the transition point 
(figure 125k ). However, at and below the transition point the behaviour is complex and 
bimodal distributions are obtained at varying locations in the layer (figure [2Sfc) ■ Similar 
behaviour is obtained for all the cases studied. The c y distributi ons bel o w the transition 
point are qualitatively similar to the distributions obtained by iMuethl (|2003) for dense 
Couette flow. 

In the region below the transition depth (y c ), the decrease in shear rate with depth, 
accompanied by an increase in shear stress with depth, implies a sharp increase in the 
apparent viscosity, defined as 7/ = Tyx/'j- Figure I2tj1 shows the variation of the apparent 
viscosity (77) with the r.m.s. velocity (it). The granular temperature is T = u 2 , thus 
figureEfilmav be considered to be analogous to the variation of viscosity with temperature 
in molecular fluids. The viscosity varies by more than three orders of magnitude and all 
the data show qualitatively similar behaviour. Near the free surface the viscosity is low 
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Figure 25. Velocity distributions at different depths the layer for 2 mm SS balls rotated at 3 
rpm. Thick solid line with filled circles represents the data at the transition point. (a),(b) Region 
at and above the transition point. (c),(d) Region at and below the transition point. 



and increases at constant u, essentially due to the increase in the volume fraction of 
solids. Below the transition point there is a gradual increase in viscosity with depth 
corresponding to a decreasing u. The variation of r\ with u in the region just below the 
transition in all cases is power law with an exponent of about —1.5. At lower values of u 
the curves turn upward and rate of increase of 77 appears to become higher. 

Figure 1771 shows the variation of the r.m.s. velocity (u) with shear rate (7). There is 
a power law dependence between the two in the region below the transition correspond- 
ing the high density region. The exponent (a) is independent of the particle diameter 
(figure 127k ). but appears to decrease systematically with increasing rpm, ranging from 
a — 1 for low rotational speeds to a — 0.5 for 9 rpm. However, data over a wider range 
are required to verify these results. Power law variation between the r.m. s. velocity and 
the shear rate have been obtained previously bv lBocauet et all (|200lj) and lMuetbl (^2003^ 
for dense Couette fl ows. The value of the exponent reported by them is around a = 0.5. 
lAzanza et all l|l999l) found a = 1 for rapid flows. 

The above results indicate the presence of two regimes in the flowing layer with a 
sharp transition point. In the upper region near the free surface, the behaviour is fluid- 
like and the velocity distributions are Maxwellian. Below the transition point the material 
appears to be an amorphous soft solid, increasing in strength with depth in the layer. 
The transition to this solid-like regime occurs at a relatively large r.m.s. velocity (~ 0.1 
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Figure 26. Apparent viscosity r\ versus r.m.s. velocity u for the different cases studied, (a) 
Open symbols represent experimental data for 2 mm SS balls. Filled circles represent data for 
2 mm brass particles rotated at 3 rpm. (b) Open symbols represent experimental data for the 
SS balls rotated at 3 rpm. Filled triangles represent data for 1 mm SS balls rotated at 2 rpm. 
Dotted lines represent a least squares fit (slope ~ —1.5) to the data in the linear region. 




Y(s"') Y(s"') 

Figure 27. Variation of the r.m.s. velocity (u) with the shear rate (7) in the the flowing layer. 
Dotted lines show least squares fits to the linear portions of the graphs, (a) Particles of different 
diameter and a rotational speed of 3 rpm. Filled symbols denote data for 2 mm brass particles, 
(b) 2 mm particles at different rotational speeds. 

m/s). Qualitatively similar behaviour is obtained for different particles and operating 
speeds. 

We conjecture that the sharp transition occurs because of the formation of a percolated 
network of particles in extended contact with each other. This is in contrast to the fluid- 
like regime where the particles interact through collisions. The contact network coexists 
with fluid-like domains and the fraction of particles which are part of the network increase 
with depth. This picture for the region below the tr ansition point is broadly consistent 
with recent measurements of lBonamv et al\ l)200 26) in which the flowing clusters were 
identified, as well as w ith non-local models based on the coexi stence of particle chains 
and fluid-like material (jMills et oi.lll999HBonamv & Miilsll2003(L 
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11 
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16 


2 
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11 


16 
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9 


14 


15 


3 


3 


7 


20 



Table 3. Measured thicknesses of the fluid-like (5f) and solid-like (S s ) layers for different cases. 



Table gives the thickness of the fluid-like (Sf) and the solid-like (S s ) regions obtained 
from the data in figures l2*2l and 12*51 The data indicates that the thickness of the fluid-like 
layer is independent of the particle diameter but depends on the rotational speed (i.e., the 
local flow rate, q). In contrast, the thickness of the solid- like region is independent of rpm 
but proportional to the particle diameter. The latter result is consistent with the finding 
that the characteristic length for exponential decay of the velocity (A) is proportional to 
diameter but independent of rotational speed. 

3.4. Comparison to theory 

IBonamv fc Millsl l|2f)0^ proposed a non-local model fo r the surface flow of granular ma- 
terial, based on an earlier work by iMills et~al\ l)l999l) . The model assumes the flowing 
granular medium as a network of transient particle chains, immersed in an assembly of 
particles behaving as a viscous Newtonian fluid. The stresses are then propagated through 
these chains in a non-local manner. The total stress is expressed as a linear combination 
of Coulombic frictional stresses, viscous stresses and the stresses due to the embedded 
chains. 

The model equations, simplified for the case of surface flows in rotating cylinders for 
low rotational speeds (oj), when inerti a! stresses can be negle cted, yield an analytical 
expression for the velocity (v x ) profile l|Bonamv fc Millsl |2P03)) 

v x (y) * - (^) (y 2 + i? 2 ) ye (-R, -S) (3.7a) 

* 7 (y + S + - (^) (y 2 + R 2 ) ye (-S, 0) (3.76) 

where A is a function of the characteristic correlation length of the transient chains. The 
shear rate (7) and the characteristic length (A) are parameters of the model. Using the 
mass balance equation at the centre of the cylinder, gives the following expression for 
the layer thickness 

Equations l|3.7|l and l|3.8|) are used to calculate the velocity profiles across the flowing 
layer. The shear rate in each case is taken to be the value obtained by fitting a straight line 
to the linear part of the experimental profile (ja, figure IT5|l . The values of the correlation 
length are found to be related to the particle diameter as X/d — 4 for all the cases. 

The predicted and the experimental profiles are shown in figure l2"51 The model predicts 
the exponential and the linear part of the velocity very well for all rotational speeds and 
particle sizes, but does not predict the flattening of the profile near the free surface. 
For the case of of 3 mm balls, the layer thickness calculated from (|3.8(l overpredicts the 
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Figure 28. Mean velocity (v x ) profi les in the flowing layer. Solid lines show the predictions of 
the model of IBonamv fc Mills ( 2003) . (a) Symbols denote data for 2 mm SS balls, (b) Symbols 
denote data for SS balls rotated at 3 rpm. 



experimental layer thickness resulting in the shift of the predicted profile (dashed line 
in figure EH 3 )- The reason for this mismatch is the experimental difficulty in correctly 
identifying the location of the free surface. If, for this case we use the layer thickness 
v alue obtained from the exper imental velocity profile, a good fit is obtained. 

lAranson fc Tsimrind l)2002(l proposed an order parameter description for the fluidiza- 
tion transition in granular flows. The order parameter specifies the fraction of the solid 
and fluid parts of the stress tensor, and is related to the number of particles in contact 
with a given particle on the average (coordination number). The dynamics of the or- 
der parameter (p) are assumed to be governed by the Ginzburg-Landau equation which 
reduces in the case of a steady, fully-developed unidirectional flow to 

where D is the characteristic diffusion coefficient for p. F(p, fj.) is the derivative of the 
potential energy density and is given by 

F(p, //) = (p - l)(p 2 - 2p*p + pi exp[-A( Ai 2 - £)]) (3.10) 

where \i = T yx /T yy and the quantities with asterisks are model parameters. IVolfson. Tsimring fc Aransonl 
( 2003) found p» = 0.6, A — 25 and /x* = 0.26 from two-dimensional discrete element sim- 
ulations for D = 2. The fluid (t^ x ) part of the total shear stress (r yx ) was obtained from 
simulations in terms of the order parameter as 

T.l x = {l-p?-\ yx = i 1 d ^ (3.11) 

where r\ is the viscosity. Simulations for 2d systems indicate r\ = 12. In the above equa- 
tions (|3.9I - 13 . 1 lfl , all quantities are made dimensionless using appropriate combinations 
of to, g and d. 

We obtain the profile of the order parameter in the layer by numerical solution of i|3.9[l 
using a shooting method with the experimentally obtained profile for the local friction 
coefficient (n(y)) as an input. The free surface is assumed to be a no flux boundary for 
the order parameter (Ddp/dy — at y = 0) and a pure solid (p — 1) is assumed far 
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Figur e 29. Comparison of the shear rate profiles predicted by the model of lAranson &: Tsimrind 
(2002) (dashed lines) to the measured shear profiles, (a) Symbols denote data for 2 mm SS balls, 
(b) Symbols denote data for SS balls rotated at 3 rpm. 

from the free surface (y — 5 a >> 1)- The profiles are independent of the domain depth 
(5a) when it is large enough. The shear rate (7 = dv x /dy) profile is then calculated from 
l|3.11b ,) using the experimentally obtained shear stress profile. 

The computed shear rate profiles are compared to experimentally measured profiles 
in figure 1291 This is a more stringent test as compared to matching the model to the 
measured velocity profiles. Predictions of the shear rate profile using the parameters 
from 2d simulations give qualitative agreement with experimental profiles. One such 
calculation for 2 mm particles rotated at 3 rpm is shown in figure 1291 (long dashed line) 
using the values given above for A, D and 77, but taking /1* = 0.4. The latter was necessary 
since no solution satisfying the boundary conditions was obtained for fi* = 0.26 for the 
particular experimental profile for /j,(y) used. We found excellent quantitative agreement 
between the model and experimental results using A = 5 and D = 2 for four of the 
cases and A = 10 and D = 2 for one case, as shown in the figure [23 However, different 
values of fi* and r\ are required to get a good fit and these are given in tabled] We found 
that n* essentially determines the depth of the layer and thus fi* was adjusted to match 
the experimental and theoretical depths. The viscosity (77) is a proportionality factor 
and we calculated it as the ratio of the maximum computed shear rate to the maximum 
experimental shear rate for each profile. Once A and D are fixed, the calculated values of 
and r\ are thus uniquely determined. The parameter A did not significantly affect the 
shape of the profiles when varied over a significant range [A £ (5, 25)), but determined the 
range of n* for which solutions satisfying the boundary conditions could be obtained. 
Similarly, changing D over narrow ranges close to D = 2 did not change the profiles 
significantly, but strongly affected the range of fx* for which solutions could be obtained. 
The graph of the order parameter variation with depth for all the cases is shown in 
figure H3 Note that according to the model, the material is not a pure fluid (p > 0) at 
the free surface (y = 0). 

The values of the parameter /i* correlate well with the effective friction coefficient (//j), 
also given in table 0] Thus fi* depends on the frictional characteristics of the material 
as may be expected. The viscosity is nearly constant for different particle sizes and 
increases with rotational speed. The latter may be expected given the increase in the 
granular temperature with rotational speed. 
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Figure 30. Variation of the order paramter (p) with depth in the flowing layer corresponding 

to the cases described in table |1] 



Material d 
(mm) 



(rpm) 



.4 



SS 



5 0.36 0.97 0.36 

5 0.40 0.85 0.34 

5 0.57 1.74 0.52 

5 0.81 2.77 0.77 

10 0.41 0.84 0.37 



Table 4. Values of fitting constants for the model of lAranson fc Tsimrind 12003) . The values 
of the effective coefficient of friction at the base of the layer, iu = fi(y s ), are also given. 



Finally, we compare the experimental results obtained to the rheological model of 
Pouliquen ljl.2|l . A direct comparison is possible in this case since /i and I can be estimated 
from experimental data. Figure ISTk shows the variation of the local coefficient of friction 
(/x) with the scaled shear rate (J) for the different data sets. The data sets are truncated 
and only the points in the dense flow region below the shear rate maximum are shown. 
The data for the lower rotational speeds are all clustered toegther while the values of 
/j, for the higher rotational speeds are significantly higher. This follows from the higher 
values of \ii for the higher rpm data seen in figure El We take fi s to the value of 
extrapolated to / = 0. The range of [i s for the low rpm data corresponds to angles 
ranging from 17 to 20 dcg., which reflec ts the error in es timation of fi. The values are in 
reasonable agreement with the results of l.Ton e.t ra/J ^OOM. Figurel3"Tb shows the variation 
of (fj, — fi s ) with the scaled shear rate (/) for 2 mm particles and for different cyindcr 
rotational speeds (3, 6 and 9 rpm). The collapse of data is reasonable, and the solid line 
is a fit to all the se ts combined. The fit yields A/i = (/j,2 — /i s ) = 0.29 and Iq — 0.49. The 
values reported bv lJop et all l)2005(l are A[i — 0.261 and Iq = 0.279, and a plot of the 
model using these values is shown in figure l5Tp (dashed line). The difference between the 
two is not very large considering the the different materials used here. Fits of the model 
gave slightly different model parameter values for 1 mm particles (A// = 0.18, Iq = 0.48) 




Figure 31. (a) Variation of the effective friction coefficient (fi) with scaled shear rate (7) for 
different particle sizes and cylinder rotational speeds as indicated in the legend. Filled symbols 
correspond to 2 mm brass balls and a rotational speed of 3 rpm. (b) Increase in coefficient of 
friction (A/i) with scaled shear rate (7) for 2 mm particles. The solid line is a fit of 1 1 .2t and 
the dashed line gives predictions of the same equation but using the parameters obtained by 
Uod el al.l <|200fltt. 



and 3 mm particles (A/i = 0.20, Iq — 0.43), but only one set of data was available for 
fitting in each case. 



4. Conclusions 

An experimental study of surface flow in quasi-2d rotating cylinders is carried out using 
streakline photography. Experiments are performed for a wide range of rotational speeds 
and particle sizes. The mean velocity profiles are linear over most of the layer depth for 
all the cases with an exponential decay near the base of the flowing layer and a flattened 
region near the free surface. The r.m.s. velocities are nearly constant in a region near 
the free surface followed by a linear decrease over part of the layer and an exponential 
decay into the bed. Further, the r.m.s. velocities show a distinct anisotropy, with the 
component in the flow direction larger by a factor of about 1.5 than that perpendicular 
to the flow direction. The number densities are nearly constant across the flowing layer 
with a sharp decrease near the free surface. The shear and the normal stresses increase 
linearly with the depth of the flowing layer initially, but the shear stress deviates from 
linearity due to wall friction. Streaming stresses are found to be significant in the region 
near the free surface. The experimental velocity profiles reported here are broadly similar 
to those in previous works, however, a wider range of system parameters is explored in 
this work. 

The experimental data is analyzed in several ways. A simple scaling is found to work 
very well. The mean vel ocity profiles for different cases scaled using the shear rate from 
lOrpe fc Khakharl <)200l[) and the particle diameter collapse to a single curve. A similar 
scaling works well for the r.m.s. velocity except for the region near the free surface. The 
shear rate obtained from a least square fit of the linear portion of the mean velocity 
profile is reasonably predicted by the model shear rate (|l.l(l but with a different value 
for the empirical constant (c). The mean velocity and r.m.s. velocity profiles near the 
base of the layer are well described by an exponential function of the form exp(— y/X), 
and the characteristic decay length (A) is proportional to the particle diameter (d) in 
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both cases, with A ~ l.ld for the mean velocity and A « 1.7d for the r.m.s. velocity. 
Thus the r.m.s. velocity decays more slowly than the mean velocity. 

The shear rate profile obtained by numerically differentiating the mean velocity profile 
data indicates that the shear rate is not constant anywhere in the layer but increases to 
a maximum and then decreases to zero with depth in the layer. A transition in the r.m.s. 
profile occurs at the same depth as the shear rate maximum: above the transition depth 
the r.m.s. velocity is nearly constant whereas below the transition depth it decreases with 
depth. Evidence of a transition is also seen in the variation of the velocity distributions 
with depth: Above the transition point the distributions are Gaussian and below the 
transition point the velocity distributions gradually approach a Poisson distribution. 

The rheology of the flow is characterized in terms of the variation of the apparent 
viscosity with r.m.s. velocity. The curves indicate a relatively sharp transition at the 
shear rate maximum, and in the region below this point the apparent viscosity varies 
as r\ ~ u . The transition observed in the measurements appears to be a percolation 
type of transition, in which there is an increasing fraction of particles below the transition 
depth which are confined, due to formation of a network of particles in extended contact 
with each other. From the viewpoint of rheology, the flow comprises two layers: an upper 
low viscosity layer with a nearly constant r.m.s. velocity and a lower layer of increasing 
viscosity with a decreasing r.m.s. velocity. The thickness of the upper layer depends on 
the local flow rate but is independent of the particle diameter. In contrast, the lower 
layer thickness depends of the particle diameter but is independent of flow rate. Thus 
after establishment of the lower layer any further increase in the local flow rate results in 
the material passing through the upper layer. The rheological characteristics of the two 
layers are different. 

Predictions of both models for t he mean velocity match the experimental velocity pro- 
files reasonably well. The model of lBonamv fc Millsl ((2003J) requires only two parameters 
: the shear rate, which can be independently obtained and the characteristic length which 
is found to be related in a simple way to the particle diameter (A = 4d). This model, 
however, does not de s cribe the low density region near the free surface. The model of 
lAranson fc Tsimrind l|2002j) describes the entire velocity profile, however, experimen- 
tal methods for independently estimating the parameter values need to be devised. A 
comp arison of the predictions of the rheological model proposed by Pouliquen l|jop et all 
2005) to experimental data also shows reasonable agreement when the increase in friction 
coeffcient with shear rate is considered. However, the coefficient of friction is found to 
be significantly higher at high local flow rates (corresponding to high rotational speeds) 
which is not in agreement with the model. The models considered show promise for 
prediction of the mean velocity profile and effective friction coefficient. More detailed 
rheological models are required for predicting r.m.s. velocities and velocity distributions 
which are important for transport. 

The results presented in this paper give a reasonably detailed picture of the rheology 
of surface granular flows in quasi- two-dimensional systems. Although the results may be 
expected to to valid in qualitative terms for 3d systems, wall effects cause significant 
quantitative differences. The walls impose two effects: one is side wall friction, which is 
accounted for approximately in this w ork, and the second is ordering b y the steric effects 
of the side walls. The measurements of lCourrech du Pont et all l|2005h indicate that wall 
friction results in the velocity at the side wall being about 20-30% lower than that at the 
middle of the gap. Our preliminary results indicate that the effects of ordering may be 
larger: the mean velocity in long cylinders is about half of that for a comparable quasi- 
two-dimensional system. Thus the ordering appears to reduce the effective viscosity of 
the system. Modelling such wall effects remains a challenge. 
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